A model of quark and lepton masses I: The neutrino sector 
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If neutrinos have masses, why are they so tiny? Are these masses of the Dirac type or of the 
Majorana type? We are already familiar with the mechanism of how to obtain a tiny Majorana 
neutrino mass by the famous see-saw mechanism. The question is: Can one build a model in 
which a tiny Dirac neutrino mass arises in a more or less "natural" way? What would be the 
phenomenological consequences of such a scenario, other than just merely reproducing the neutrino 
mass patterns for the oscillation data? In this article, a systematic and detailed analysis of a model is 
presented, with, as key components, the introduction of a family symmetry as well as a new SU(2) 
symmetry for the right-handed neutrinos. In particular, in addition to the calculations of light 
neutrino Dirac masses, interesting phenomenological implications of the model will be presented. 
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I. INTRODUCTION 



There are strong indications- the latest of which came from the SuperKamiokande collaboration Q- that neutrinos 
do have a mass, albeit a very tiny one, and, as a result, "oscillate". The exact nature of the masses as well as the 
oscillation angles is an important subject which is under intense investigation [0. Consequently, there exists many 
interesting models which, in one way or another, try to accomodate most of the known data. It is perhaps prudent 
. to think that the subject of neutrino masses and oscillation is still a very open one. 
(f-) ' It is fair to say that the extreme smallness of neutrino masses suggests something very peculiar about these particles. 

This peculiarity could come from the way the neutrinos obtain their masses and/or from the very special nature of 
■ the neutrinos themselves which distinguish them from all other particles. For example, do right-handed neutrinos 
(present in most models of neutrino masses) carry quantum numbers which are absent in some or all other (left- or 
r"| ■ right-handed) fermions? After all, right-handed neutrinos, if present, would be singlets under SU(3)®SU(2)l®U(1)y 
O ^ anyway. 

JL . Most efforts on the problem of neutrino masses, at least on the model-building front, are concentrated on the 
construction of lepton mass matrices based on various ansatzes. There is one common assumption present in many 
of such models, which is one in which light neutrino masses arise from a see-saw mechanism [||. The smallness of 
neutrino masses would come from an expression that goes like m 2 D /M, where mri is a Dirac mass , and Ai is a 
. £^ | Majorana mass which typically is very much larger than mry . In these models, the scale of new physics M, as 
suggested by the lightness of neutrino masses, would be some kind of Grand Unified scale or even the breaking scale 
of Left-Right symmetry models [0. (Lepton number is not a conserved quantity in this class of models.) The see-saw 
mechanism is a very elegant approach which is widely embraced. 

However, one could not help but wonder if there might be some other mechanism for obtaining tiny neutrino masses, 
and if so, how it would fare compared with the see-saw mechanism. Would this new mechanism shed light on other 
important issues? What would be its scale of new physics? Can one find an experimental distinction between the two 
mechanisms? This was the topic discussed in [El. 

At the present time, it is not clear that, if neutrinos do have a mass, it would be of the Majorana or Dirac type. 
As we have mentioned above, with Majorana neutrinos and the see-saw mechanism, one could "easily" obtain small 
neutrino masses. Now if the mass were to be of the Dirac type, one can straightforwardly write down a gauge- invariant 
Yukawa coupling in the SM itself (endowed with right-handed neutrinos, of course). But to obtain a small neutrino 
mass, one has to put in by hand a Yukawa coupling which is incredibly small, of the order of 10 . Such a fine 
tuning is highly unnatural and that might be the reason why little attention is given to the construction of models 
based on Dirac neutrino masses. Did we leave something out by ignoring it? What if the mass is truly of the Dirac 
type? Until this question is settled, it is worthwhile to investigate possible alternatives to the see-saw mechanism. 
This paper and a previous one || propose one of such alternatives by constructing a model of Dirac neutrino masses 
where the smallness of their values arises dynamically. One of the criteria used in building such a model is the wish 
to go beyond the mere presentation of a neutrino mass matrix. In particular, we would like to see if there might be 
other phenomenological consequences which could be testable: New particles, new physics signals, etc.. This is the 
aim we set about in building our model. 
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The construction of the model presented in || was based on the following questions: If neutrino masses were so 
small compared to all other known masses, would there be an appearance of a special symmetry when one lets the mass 
go to zero? Could this special symmetry, if it exists, be a peculiar feature of the right-handed neutral leptons alone? 
Could there be additional purposes for its existence other than providing a small mass for the neutrinos? In other 
words, can one learn something more from it? It was found in Ref. || that there is indeed an interesting symmetry 
which acts only on the right-handed neutrinos and which, in addition to providing a reason for the smallness of the 
neutrino masses, also constrains the nature (even or odd) of the number of generations. Furthermore, the way in which 
neutrino masses are constructed can be used to build a model for charged lepton and quark masses. In addition, this 
particular way of constructing masses might even have some bearing on the strong CP problem. Last but not least: 
Are there additional tests of various neutrino models other than neutrino oscillations? For the see-saw mechanism 
with Majorana neutrinos, one already sees that one of such additional signals is, for example, the phenomenon of 
ncutrinolcss double beta decay. As it will be presented below, the addtional signals of the model presented here will 
involve a number of very concrete predictions: the absence of ncutrinolcss double beta decay, the possible presence of 
"low mass" ( a couple of hundreds of GeV e.g.) vector-like fermions, among other things. In particular, the detection 
of these vector-like fermions do not in any way involve neutrinos. 

One particularly important feature of our model is the following predictions for neutrino oscillations, assuming only 
the validity of the atmostpheric and solar neutrino data: 1) The three light neutrinos are nearly degenerate; 2) If the 
light neutrinos have a mass large enough to form a component of the Hot Dark Matter (HDM) || then only the MSW 
solution to the solar neutrino oscillation is favored; 3) If the vacuum solution to the solar neutrino problem turns 
out to be the correct one, our model will only be able to accomodate tiny neutrino masses, around !Q~ 3 eV or less, 
ruling out near-degenerate neutrinos as components of HDM. As a result, in our model, one cannot have both vacuum 
solution and HDM. We will show below the correlation between the masses and the differences of mass squared, Am 2 , 
which enter the neutrino oscillation phenomena. 

Assuming the existence of the aforementioned symmetry, how can one construct Dirac neutrino masses to be 
dynamically small? By "dynamically", it is meant that the mass is zero at tree level and that any non-zero value 
would have to arise at the one-loop (or more) level. Now, the peculiar (and toughest) thing about neutrinos is the fact 
that their mass is so small- at least eleven orders of magnitude smaller than the electroweak scale. In constructing 
our model for Dirac neutrino masses, it is then reasonable to ask under what conditions would the dynamical Dirac 
mass of the neutrinos obtained at the one loop level be "naturally" small, i.e. devoid of excessive fine tuning. In this 
paper, we present the following interesting results: In the four-generation model, it is found that the fourth neutrino 
can be naturally heavy while the other three obtain their masses at one loop, with the result that these masses can 
be tiny provided some ratios of masses of particles which participate in the loop integration are "large" , regardless 
of their actual values. This is interesting because, as we shall see below, some of the particles which participate in 
the loop integration, in particular the lightest ones, can have masses as low as a few hundred GeVs and which could 
provide a direct test of this model. We will also see that, in order to obtain very small neutrino masses, at least one 
of the particles needs to be much heavier than the lightest one- a result which is somewhat reminescent of the see-saw 
mechanism. We will also see that the mass of the light neutrinos is intrinsically tied to the extra global symmetry 
present in the scalar sector of the model. In fact, the extra Nambu-Goldstone (NG) bosons which are not absorbed 
by the (family and SU(2) UR ) gauge bosons acquire a mass due to the presence of the gauge-invariant "cross-coupling" 
terms in the potential which explicitely break the extra global symmetry. 

The above brief statement will be made clearer in the discussion of neutrino masses. Notice, in particular, that the 
result given for light neutrino masses in |5| is only a very special case of the present discussion. 

The plan of the paper is as follows. First, the model is presented with a description of the gauge structure along 
with its particle content. It is shown how a new symmetry prevents neutrinos from obtaining a mass unless it is 
broken. Next, the special properties of this extra symmetry associated with the right-handed neutrinos are discussed. 
In particular, if that symmetry is a chiral SU (2) as is the case in this paper, nontrivial constraints coming from 
the nonperturbative Witten anomaly Q can be applied to the nature of the number of families. This is the extra 
bonus mentioned above. The paper then proceeds to discuss the generation of light neutrino masses, principally by 
radiative corrections of the type mentioned above. It is then followed by a discussion of the neutrino mass matrix. 
In particular, we will present the correlation between the values of the neutrino masses and Am 2 . Most importantly, 
we will show how Am 2 increases or decreases with the masses themselves, with two resulting implications: either one 
has HDM and MSW or vacuum solution and no HDM. Either of these solutions will have an important cosmological 
implication. We end the paper with a brief discussion of the charged lepton mass matrix, the primary purpose of 
which being the wish to complete the discussion by presenting some examples on what the oscillation angles might 
look like. A followed-up paper will deal seperately with the charged lepton sector and, as a consequence, with a full 
discussion of the angles. 

We would like to emphasize for the purpose of clarity that the charged lepton sector (which will be dealt with in a 
separate paper) is different in structure from the neutrino sector, as we shall see below, and does not have the same 
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hierarchical structure. The fact that, in this model, the three light neutrinos are nearly degenerate does not imply 
that it would be the same in the charged lepton sector. In fact, it is not as we will show in a subsequent paper. 

Finally, a section will be devoted to various other phenomenological implications of the model. 

We shall assume throughout this paper the existence of right-handed neutrinos. 

Since this manuscript is meant to be comprehensive, and hence lengthy, one could skip the three subsections of 
the next section, after first reading its introduction. (Its reading is nevertheless recommended because the physics 
motivations are discussed there.) 

II. A MODEL 

It is well-known that all that is needed to give neutrinos a mass is to simply add extra right-handed neutrinos to the 
Standard Model. One can then construct a (Dirac) mass term with an arbitrary Yukawa coupling, q v Il ^vr + H.c, 
which can be made to be as small as one wishes. This, of course, is unsatisfactory because, if neutrinos have masses 
in the eV range or less, this would require the Yukawa coupling, g v , to be of O(10~ n ) or less. Fine tuning to such a 
precision is normally considered to be unnatural. At this point, one might be tempted to try to explain this fact by 
simply invoking a fourth generation with a democratic mass matrix, at least for the neutrinos, as has been done by 
Ref. Q . The diagonalization of the neutrino mass matrix would then give one heavy eigenstate and three massless 
states. By adding some arbitrary phases to the mass matrix, one can "provide" a small mass ( depending on the 
values of those phases) to the three neutrinos. This purely phenomenological ansatz (Ref. ||) appears to "fit" the 
recent data on neutrino oscillations with the appropriate choices of the phases. However, the fourth generation lepton 
masses came out to be extremely heavy and split, which practically seems to be ruled out by analyses of precision 
experiments I|. 

In ||, a model of Dirac neutrino masses was constructed and based on a four generation scenario that was very 
different from the democratic ansatz made in ||. One of the reasons for using such a scenario is the fact that, as of 
the present time, a fourth generation is not ruled out by experiment and, as a consequence, it is interesting to explore 
its possible implications. A recent review [ p"l| ] gave a comprehensive discussion of various topics concerning quarks 
and leptons beyond the third generation, including the present experimental status and future searches. 

If a fourth generation were to be used in the investigation of neutrino masses, one should keep in mind various 
phenomenological constraints concerning not only leptons but also quarks. For instance, constraints on the p parameter 
limit the mass splitting within each doublet of extra quarks and leptons: the up and down members of a fourth 
generation should be very close in mass. They should be long-lived enough to escape present detection. This, in 
turns, tells us something about the mixing between the fourth generation and the other three. All of these issues have 
been described in jll]]. In the construction of the model presented in Ref. ||, these phenomenological constraints 
were kept in mind. 

As mentioned briefly in the Introduction, our approach, as described in Ref. ||, is based on a dynamical justification 
for for the small value of the neutrino Yukawa couplings. The question that was asked was: Could there be a scenario 
in which a symmetry appears as one lets the Yukawa coupling go to zero? The tiny Yukawa coupling which would give 
the neutrino a very small mass would then arise dynamically when that symmetry is broken. These Yukawa couplings 
then appear as effective couplings which could be small for dynamical reasons and are not fundamental parameters 
that are put in by hand and which are needed to be fine tuned. What is the nature of that symmetry and how a 
dynamical Yukawa coupling appears will be the subject of this section and the following two. 

It is obvious that an extension of the Standard Model(SM) is needed in addressing the above issues. One simply 
cannot stay solely within the SM if one wishes to deal with the mass of the neutrinos. What it is that one needs when 
one goes beyond the SM is a matter of taste, modulo a very obvious requirement: predictability of new phenomena 
or particles which can be tested. 

We first describe the model, presenting its gauge structure and representations. Next, explanations are provided for 
the reasons behind the choices of the extended gauge group and its particle content. The crucial assumption here is 
the existence of two new symmetries, one of which will be particular to right-handed neutrinos, as alluded to earlier, 
and the other one is a family gauge symmetry. As we shall see below, it is the breaking of these new symmetries that 
will give a mass to the neutrinos. 

In this work, the SM is extended in the following way. Generically, it takes the form: SU(3) C ® SU(2)l <8> U(1)y ® 
(Family symmetry) ® (right — handed neutrino special symmetry) . Why a "Family symmetry"? This is so for two 
reasons: a) We wish to investigate the family replication problem and the mixing among different generations; b) 
The special symmetry endowed by the right-handed neutrinos might have some bearing on the family symmetry 
itself. After all, if one would like to investigate the family problem, some kind of family symmetry-be it discrete or 
continuous, global or gauge- is needed. Why a special symmetry for the right-handed neutrinos? The reasons were 
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already expounded above: To provide a framework for an understanding of the smallness of neutrino masses. Our 
next task is then to determine what this special symmetry might be and what form the family symmetry might take. 
Our model is described by: 

SU(3) C <8> SU(2) L U(1) Y <8> SO(N f ) <g SU(2) UR (1) 

where SO(Nf) and SU(2) VR are the family gauge group and the special gauge group for the right-handed neutrinos 
respectively. The particle content of the model is listed in Table 1. Notice that we have denoted the right-handed 
neutrinos by vr = (vr^r) because they are assumed to transform as doublets under SU(2) VR . The two options 
listed for the right-handed neutrinos as well as the meaning of the non-standard particles will be discussed below. 
We would first like to explain the choices of the extra gauge groups. Here, the extra symmetries are chosen to be 
gauge symmetries because, as it is well known, powerful constraints can be obtained from models built on the gauge 
assumption. 

A. Why SU(2) VR 1 

Let us first look at Table 1. In this model, all standard (left-handed and right-handed) particles are singlets under 
SU(2)y Hence the subscript vr. In this respect, SU(2) VR is very different from SU(2)r of the popular Left-Right 
model 11 . In that model, right-handed quarks and leptons form doublets under SU(2)r, for every family. Because of 
our assignment, all weak interactions among standard particles are pure V-A, in contrast with the Left-Right model. 
What is the motivation behind our choice that makes it so different from the Left-Right model? To answer that 
question, let us recall an interesting feature of chiral SU(2): the presence or absence of the so-called Witten global 
anomaly. 

If chiral fermions transform as doublets under SU(2), there exists a nonperturbative anomaly- the so-called Witten 
anomaly 0- associated with an odd number of doublets. Briefly speaking, this is so because the fermionic determinant 
y/deti ft (A,j,) changes sign under a "large" gauge transformation A 1 ^ = U~ 1 A fl U — iU~ l dfJJ if the number of 
chiral doublets is odd. This would make the partition function Z vanish and the theory would be ill-defined. This 
nonperturbative anomaly would then require the number of Weyl doublets to be even in order for the theory to be 
consistent. (This ambiguity in sign stems from the fact that the fourth homotopy group TLl(SU(2)) — Z%.) Other 
groups that also have similar non-trivial constraints are Sp(N) for any N and O(N) for N < 5. 

It is amusing to recall a well-known but forgotten fact about the SM. There the chiral gauge group is SU(2)l. Each 
family contains one lepton and three quark doublets and, as such, is free from the global Witten anomaly. (Let us 
recall that the cancellation of the perturbative triangle anomaly in the SM only relates the lepton charge to that of 
the quark.) If, instead of three, the number of colors, N c , were arbitrary, the freedom from such an anomaly would 
require 1 + N c to be even, and hence, N c to be odd, namely N c — 3, 5,... Why nature choses N c = 3 instead of some 
other odd number is a question which can only be answered in the context of some deeper theory such as e.g. SU(5). 
Although the Witten anomaly does not fix the size of N c , it is nevertheless a powerful constraint in the sense that, 
once a fermion content is known (e.g. one color singlet (leptons) and one fundamental representation (quarks) in the 
SM), N c is constrained (e.g. odd in the case of the SM). 

The above simple lesson taught us something about the powerful constraint that a chiral SU (2) exerts on the 
number of chiral doublets. This is the reason why it is chosen to be the special symmetry of the right-handed 
neutrinos. Let us contrast the constraint coming from SU(2) UR with that coming from SU(2)r (Left-Right model). 
For our model, with SU(2) VR1 only tjr transforms as doublets. Absence from the Witten anomaly then requires the 
number of such doublets to be even. If t]r carries, in addition, family indices then the anomaly requirement restricts 
the number of generations to be even such as in Option 1 as indicated in Table 1. If there exists an vr which is a 
family singlet (denoted by t]'r), the number of generations would be odd such as in Option 2 of Table 1. With the 
Left-Right model, each family contains four doublets of SU(2)r\ (vrtCr) and (ur^r)i with i = 1,..,3. Therefore, 
the Witten anomaly requirement is automatically satisfied per family. This is one of the few differences between our 
model and the Left-Right model. 

A final word of caution is in order here. Although the Witten anomaly constraint allows us to make a statement on 
the evenness or oddncss of the number of generations- a subject to which we shall come back in the next subsection, 
it does not determine that number. This should come from a deeper and as-yet-unknown theory. Our goal is much 
more modest: Given a fermion content (Option 1 or 2 below), we can say whether or not the number of generations 
is odd or even, and that is all. We shall however try to constraint that number from a different route which is more 
phenomenological in nature, and point out the differences between Option 1 and 2. 
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B. Why SO(Nf)? 



In the construction of any model, there is a time- honored requirement: the absence of the perturbative triangle 
anomaly. Even if the Witten anomaly were absent, this requirement is a must for any gauge theory. (It just happens 
that, in the SM, both requirements are simultaneously satisfied.) In our case, if a family index is assigned to all 
standard fermions and to r]n, the family gauge group that is chosen cannot be a vector- like theory, which is anomaly- 
free, because rju posseses an additional quantum number, that of SU(2) VR . This is unlike QCD or even the Left-Right 
model if left and right-handed fermions carry similar family quantum numbers. A safe group and representations have 
to be chosen. 

The choice made in this paper is SO(Nf) for the family gauge group, with chiral (left- and right-handed) fermions 
transforming as (real) vector representations with Nf components each. As such, the model is also free of the 
perturbative triangle anomaly. 

Our model based on SU(3) c ®SU(2) L ®U(l) Y ®SO(N f )®SU(2)„ R with an even number of SU(2) VR - doublets and 
chiral fermions transforming as vector representations of SO(Nf) is free from both nonperturbative and perturbative 
anomalies. 



C. Constraints on Nf 

As shown in Table 1, there are two options for tjr, each of which should contain an even number of SU{2) 1Jll 
doublets. 

a) Option 1: 

rj R carries the family index a — l..Nf where Nt = 2, 4, 6, 8, ... 

b) Option 2: 

Here we have rj' R (a family singlet) and r/ R . The constraint is now 1 + Nf = even, which means that Nf = 3, 5, 7, .. 
(excluding the trivial case of 1 family). 

Unlike the SM where one knows the fermion content for each family, i.e. quarks and leptons, and hence the nature 
of N c - it is odd- our scenario involves incomplete experimental informations, and as such, the nature (odd or even) 
of Nf cannot be completely fixed. Each choice, however, represents a distinct particle content (no family singlets for 
the even option and one family singlet for the odd option) which implies a possible distinct route for a yet-unknown 
unification. 

Recognizing the fact that there are deep differences between the even and odd options-a point to be discussed below- 
and in the absence of a deeper theory, one might wonder what can be done to narrow down the choices, not between 
odd or even, but within each option itself. Below we present an argument that could help in finding a way to further 
restrict Nt. This argument is only suggestive, being a combination of "theoretical prejudice" and phenomenological 
constraint. 

One might require that gauge couplings are free from Landau singularities below the Planck scale in such a way that 
unification of the SM gauge couplings, if it exists, occurs in the perturbative regime JTo[ ] . With this criterion, one can 
see that the even option can only accomodate Nf — 2, 4, 6, while the odd option can only accomodate Nf = 3, 5. This 
is because for Nf > 7, one or more gauge couplings will "blow up" before the Planck scale. There are no reasons, in 
the absence of a deeper theory, to rule out any of the above choices. This will require other yet-unknown conditions. 
The only thing one can say, in the context of our model, is that electroweak precision experiments appear to rule out 
Nf > 5 U and and that existential facts tell us that Nf is at least three. This leaves us with the choice Nf — 4 for 
the even option and Nf = 3 for the odd option. 

If Nf < 4 comes from the above argument, what then is the role of the Witten anomaly in all of this? It tells 
us about the particle content of the right-handed neutrinos. For Nf = 4, the right-handed neutrinos are simply 
rju = (1, 1, 0, 4, 2) while for Nf — 3, one has rj R — (1, 1, 0, 3, 2) plus a family singlet r)' R = (1, 1, 0, 1, 2). What observed 
differences can there be between these two options? The former predicts the existence of a fourth generation whose 
consequences have been recently discussed in Refs. |l(]] and |ll],|l2|]. The latter predicts the existence of a neutral 
family-singlet ij R (doublet under SU(2) UR ) which could have cosmological consequences of a yet-unknown nature. In 
addition, as we point out below, it appears that the even option prefers three almost degenerate light neutrinos while 
the odd option prefers a hierarchical structure for the light neutrinos. If a fourth generation is discovered, which alone 
does not necessarily imply the even option presented here, and if the light neutrino masses are convincingly "proven" 
to be nearly degenerate (instead of a hierarchical structure), the even option might be viable. Furthermore, as we 
shall see below, another possibility for testing this model is to look for signals of some of the lightest particles- the 
vector-like fermions- which participate in the loop diagram of Fig. 1. As discussed below, the light neutrino masses 



5 



depend only on the ratios of these masses and not on their magnitudes and these vector-like fermions can be as light 
as a few hundred GeVs. 



III. NEUTRINO MASSES 



This section will be devoted to the discussion of how neutrino masses can be generated in our model for Option 1 . 
We shall comment on Option 2 at the end of the manuscript. We shall concentrate only on the lepton sector and, in 
particular, on the neutrino one, leaving the full discussion of the charged lepton and quark sectors for a subsequent 
publication. 

Since we will be dealing only with Dirac neutrino masses, we shall require that all fermions be endowed with a 
global B — L symmetry. Since we are concerned only with leptons in this section, a global L symmetry is sufficient 
for the present purpose. This global L symmetry would prevent a Majorana mass term of the type rf R a rii aR , where 
i = l,2 and a = 1, ..,4. Only Dirac masses will be allowed. 

There might be other suggestive reasons as to why Dirac masses for the neutrinos might be attractive. For example, 
a combined fit of massive neutrinos as components of Hot Dark Matter (HDM) and atmospheric neutrino oscillations 
seems to prefer a scenario in which two or three light neutrinos are nearl y d egenerate and have mass in the O(eV) 
range. Recent data on neutrinoless double beta decay (or absence thereof) appear to rule out Majorana neutrinos 
heavier than 0.2 eV, at least in the simplest versions. Here it will be shown how, in our scenario, one can obtain three 
near-degenerate neutrinos whose mass can be of the order of a few eV's and is of the Dirac type. Consequently, in 
our model, there will be no neutrinoless double beta decay, and hence no contraint on the Dirac neutrino masses from 
such a search. 

As we have discussed above, Option 1 contains no family-singlet fcrmion field and freedom from the Witten anomaly 
dictates that the number of families should be even. Furthermore, we have argued that this even number should be 
four. As a result, the gauge group for this option is: 

SU(3) C <g> SU(2) L ® U(1)y ® SO(A) ® SU{2) VR (2) 

The reader is referred to Table 1 for a list of particles that participate in this model. 



A. Computation of the diagonal elements of the 4x4 neutrino mass matrix 

Without the extra vector-like fermions, F, Mi and M 2 , the only gauge-invariant Yukawa coupling involving leptons 
would be Cy = gEl < L<t>e a R + H.c, (where a = 1, ..,4 is the family index), giving rise to equal masses for the charged 
leptons. Unbroken SU(2) u r forbids a similar term for the neutrinos and they remain massless at this level. (Notice 
that, since we are only interested in Dirac neutrino masses, a gauge- invariant Majorana mass term of the type rf R a r\i a r 
is forbidden by L symmetry.) We know that the charged leptons are not degenerate in mass. We also know that 
the width of the Z boson |L4j constrains the mass of the fourth neutrino to be larger than half the Z mass. This is 
where the vector-like fermions listed in Table 1 come in. Because of their vector-like nature, they can have arbitrary 
gauge-invariant bare masses. It is seen below that some of these masses can be as low as a few hundreds GeVs and 
are thus accessible to future experimental searches. 

The Yukawa part of the Lagrangian involving leptons can be written as 

^lepton = 9El a L 4>e aR + Gj a L n a F R + G Ml F L <j>M 1R + Gm 2 F l 4>M 2 r + G 2 Mi L Cl a e R + 

G 3 M 2L pXr + M F F L F R + M 1 MilM 1r + M 2 M 2L M 2R + h.c. (3) 

The assumption of an unbroken L symmetry forbids the presence of Majorana mass terms as mentioned above. 

Notice that the values of Mp,i, 2 are arbitrary. What they might be will be the subject of the discussion presented 
below. After integrating out the F, A4i, and M. 2 fields, the relevant part of the effective Lagrangian below Mp,i,2 
reads 

1%SL = 9El^e aR + G E Tl{£l a <t^)e pa + 

G N l2(nJp?)4 R + H.c, (4) 

where 

GE ~ M F Ah ' GN ~ M F M 2 ■ (5) 
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This is a tree-level effective Lagrangian whose consequences are now presented. 

Let us discuss the implication of each term on the right-hand side of Eq. (^). As stated in the preceding paragraph, 
the first term gives rise to equal masses for the charged leptons. The second term would lift the degeneracy of the 
charged lepton sector once f2 acquires a vacuum expectation value (VEV). The third term gives rise to a neutrino 
mass once both fi and p acquire a VEV. It is clear that, in our model, neutrino masses can appear only when both 
50(4) and SU{2) UR are spontaneously broken while the charged lepton masses are non zero (but equal) even if SO (4) 
is unbroken. Only when 50(4) is broken will the charged lepton mass degeneracy be lifted. 

Let us assume: < Ci >= (0, 0, 0, V) and < p >= (0, 0, 0, V <g> si), where si = ( J ^ . Notice that each component 

(under 50(4)) of p transforms as a doublet under SU{2) VR . If we denote the 4th element of t]r by (Nr, Nr), one 
can use the above two VEV's along with < <f> >= (0, v/V2) ( v ~ 246 GeV) in Eq.(|4|) to write down a Dirac mass 
term for the 4th generation neutrino, namely 

G N -^=N L N R + h.c; G N = GxG M2 G 3 ^ V , (6) 

V2 Nip M-x 

giving 

V 

m N = G N — . (7) 

At tree level, all other neutrinos are massless. Their masses arise at the one-loop level as shown below. The Dirac 
mass of the fourth neutrino could be rather heavy. In fact, it is not unreasonable to expect Gx, Gm 2 an< i O3 to be of 
the order of unity. In consequence, as long as 

VV'/MfM*s ~ O(l), (8) 

one might expect the fourth neutrino to be even as heavy as 175 GeV. Certainly, the LEP bound of Mz/2 can easily 
be satisfied. 

Why are the other three neutrinos massless at tree level? Firstly, it is so because, from Eq. (S) and Eq. (|J), one 
can see that, after integrating out the heavy vector-like fermions, there is no (tree-level, dimension 6) operator which 
contains, as a factor, a term such as l™4>ViR' wnere m = 1,2, 3 is a family index, which would give rise to a mass 
term for the three light neutrinos. An effective (dimension 6) operator which contains the aforementioned term would 
necessarily come from a loop integration such as the one shown in Fig. 1. Just like the various terms which appear 
in Eq. ([!]), this effective operator would also contain the scalar fields f2 and p. It would appear as 

iT4>vr R (n a p ai ). (9) 

As pointed out in the Appendix, a term such as (fl a p a i) appears as part of a quartic term in the potential which 
explicitely breaks the extra global symmetry that the scalar sector posesses. As a result, the extra NG bosons are, 



in fact, pseudo NG bosons and acquire a mass which is proportional to the coupling A4 as shown in Eq. (A8) of the 
Appendix. 

In order to compute the one-loop contributions to neutrino masses, let us recall, in this section, the results obtained 
in the Appendix concerning the relevant mass eigenstates in the scalar sector. We have 

i?4 = cosai?4 — sin a/14, (10a) 
hi = sin aH^ + cos a/14, (10b) 
Oj = cos/3f2i — sin (3Repi, (10c) 

Re pi = sin/Jfij + cos (3 Repi, (10d) 

where i = 1,2,3 and where the states with the ~ sign are mass eigenstates. The Yukawa couplings which will be 
involved in the computation of neutrino masses can now be written in terms of the mass eigenstates. For example, 
GilfVl a FR can be written as 
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dlln^FR = GiZ|(cos aH 4 - smah 4 )F R , (11) 

Ch? L SliFR = GiZl(cos/3f2i - sux0Repi)F R . (12) 
where z = 1,2, 3. Also, GsM^Lp^if^R ( TO — 1) 2) can be now written as 

G3M2LP1V4R = G 3 M 2 L(sin aiJ 4 + cos a/14 + iImp 4 )\r\\ R , (13) 

G 3 M 2L p\ v l R = G 3 M 2L (sin + cosfiRepi + ilmp^r^ . (14) 

The above equations, in addition to G]tf 2 Fi,(f>M-z R , form the basis for constructing the one-loop diagrams as shown in 
Fig.l. As one can immediately see, the only scalars that participate in the loop integration are Hi, hi, fij, and pi. 
The contributions to the light neutrino masses will contain a factor cos/? sin j3 — sin(2/3)/2 for fli and -cos (3 sin (3 for 
Repi. 

The masses of the p hysical Higgs scalars, H4 and hi, and those of the pseudo NG bosons, Repi (i — 1,2,3), are 
given by Eqs. ( A4| A8| ) in the Appendix. Since the one- loop contributions to the 4th neutrino mass are expected to 



be small compared with its tree-level value, we shall concentrate in this section on the light neutrino masses. There 
we shall be concerned only with fii (NG bosons) and Repi (pseudo NG bosons) (i = 1, 2, 3). In the 'tHooft-Feynman 
gauge, the NG bosons will have a propagator with a mass which is that of the family gauge bosons. We shall denote 
it by Mg- We shall call the mass of the pseudo NG bosons, Mp. 

The result obtained from the diagrams as shown in Fig. 1 for the three light neutrinos is 

v 

m v = G v — , (15) 
v2 



where 



and where 



G v = G 1 G M2 G 3 ^^-{I{n)-I{Rep)), (16) 

32 7T Z 



1 M F [MUM G H^)-M P H^)) + M G M P H^)] 
^ ~ ^ = Mr-M^ (M G -M^M p -M F ) °- - ^ - M ^ ^ 

For notational convenience, we shall define: 

AI(G,P) = I(Q)-I(Rep), (18) 

It is convenient to express the mass of the light neutrinos by the following ratio: 

m u M F M 2 sin(2/3) 



m N VV 32 7T 2 



■AI(G,P), (19) 



where tojv is defined by Eq. (^). 

One should mention for completeness the tiny one-loop contribution to the 4th neutrino mass. If we denote by this 
contribution by 61714, it is straigthforward to see that it is given precisely by the same formula for the light neutrino 
mass, Eq. (|T5|), with the following replacements: — > a, Mq —> M Ri , Mp — > M^ 4 , namely 

8UI4 = Gi^=, (20a) 
v 2 



G4 — G1GA/2G; 



sin(2a) 

32 7T 2 



AI(G,P), 



(20b) 



where the form of /(i?4) — I(hi) is identical to Eq. (17) with the replacements as mentioned above. This contribution 
will play an insignificant role in the mass matrix, but it has to be mentioned for completeness. 
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The above results were obtained at one loop. One wonders if higher loop contributions might be significant. It 
turns out that, because of the nature of the interactions, the next correction occurs at the three loop level. It means 
that the correction to the one-loop light neutrino mass is at the two-loop order. Considering that already the one-loop 
result is 0(< 10~ 10 ), a two-loop correction to that result would most likely be insignificant, even for the mass splitting 
to be discussed below. Above all, the experimental results are far from being precise enough to even contemplate such 
a tiny correction. From hereon, we shall assume that these three-loop corrections are insignificant in the computation 
of the mass splittings. 

At this stage, the three light neutrinos are degenerate. A discussion of the lifting of the degeneracy will follow a 
more general discussion of the implications of Eq. (|l9|) . It is clear that the "light family" symmetry would have to be 
broken in order for the "light" fermions to mix. It is also clear that the neutrino masses (one heavy and three light) 
derived so far represent only the diagonal elements of a 4 x 4 neutrino mass matrix. If the discussion presented in this 
section on light neutrino masses is to be at all interesting, it is imperative to assume that the bulk of at least one, if 



not all, of the light neutrino masses comes from Eq. (15) 



At this point, an important remark is in order here. As we have stressed above, the near-degeneracy of the light 
neutrinos in no way implies that a similar situation will occur in the charged lepton sector. In fact, we will show in a 
separate paper that this will not be the case. 

Under what conditions will G v be of the order of 10 -11 or less? First of all, as we have seen from Eq. (0), in order 
to have a "heavy" fourth neutrino, one should have G1GM2G3 m^m^ ~ ^(-0- This puts a condition on the angle (3 
itself, namely (tan/? = V'/V) 

tan/3~ 1 M F M 2 
G1GM2G3 V 2 

As we have stated earlier, it is not unreasonable to assume that G\, Gm 2 an d G3 to be of the order of unity. With 
Mq ~ g 2 V 2 (where g is the SO (A) gauge coupling), Eq. ( ^l|) becomes 

9 M F M 2 , s 

^~ 9 W G W G (22) 

The above estimate for the constraint on the angle (3 will be used in our computation of the light neutrino masses. 
With this in mind, we can now proceed to make an estimate of the ratio m y /mjv, where now V V /Mp M 2 ~ 0(1) 



and Eq. ( |19| ) becomes 

m„ sin(2/3) 



771 jy 32 7T 
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(7(A) - I(Rep)). (23) 



As we have seen above, the result (|23J) depends only on ratios of masses of the particles in the loop integral and not 
on their absolute values. Because of that fact, the results will be shown in units of Mp which can be as small or as 
large as one wishes. 

Before moving on to discuss the implications of Eqs. ( |l9|) and (p3|), one remark is in order here. From Eq. jl^), 
one can see that the light neutrino mass vanishes when Mg = Mp. Since there is no reason (as far as the present 
construction of the model is concerned) for this equality to be valid, we shall dismiss this possibility. We shall 
concentrate instead on the criteria for having small m v for arbitrary Mq and Mp (and Mp and M 2 as well) . 

The results are shown in Figs. 2, 3, 4 and 5. A few comments are in order here. First of all, as we have mentioned 
above, our results depend on ratios of the four masses which enter the loop integral: Mp, Mp, Mq, and M%. One can 
symbolically denote one of the masses as M = 1, and the other three will be multiples of that chosen one. Which one 
should be chosen is a matter of convenience and phenomenological interest. In particular, we choose Mp = 1 because 
there is a possibility that the vector-like fermions F could be detected if their masses are low enough. 

A glance at Figs. 2-5 reveal that it is relatively easy to obtain a very small ratio R = m v /mj^. In particular, one 
can see that large values of M 2 , the mass of the singlet fermion field M. 2 , are sufficient to obtain small values for 
R = m v /mN. For instance, one can see that, roughly speaking, R = m^/m^ < 10~ n when M 2 > 10 6 (in units of 
Mp). Although conceptually quite different, the above fact is very reminescent of the see-saw mechanism in that there 
is one large scale: Majorana for see-saw, M 2 for this scenario, and one "small" scale: Dirac mass mp, for see-saw, Mp 
for this scenario. The important point that we wish to make is the fact that the general result obtained here, namely 
the smallness of light neutrino masses, does not depend on one particular combination of masses which would imply 
fine tuning, a point which was not made quite clear in Ref. Q, but only on "large" ratio of masses whatever they 
might be. In this sense, the smallness of neutrino masses in our scenario is no less natural than the ones obtained 
from the sec-saw mechanism. 







In Figs. 2-5, we show the results for the case M 2 > Mq- There is, of course, absolutely no reason for this ordering. 
It is a matter of presentation. We obtain exactly the same results with the roles of M 2 and Mq reversed. As 
can be inferred from the figures, for a given value of Mp (Mp = 1), R = m v jm^ < 10~ n if the ratio M0/M2 
is below a certain value. For example, for Mq < 10 5 , one has Mq/M 2 < 10~ 3 , while for Mq > 10 7 , one has 
Mq/M 2 w 10~ 2 — 10 _1 . What this says is that the larger the mass is (e.g. Mg), the less mass splitting is needed in 
order to have a small R. 

At this stage, we can only say that m v can be very small. What we cannot say is exactly what its value should be. 
This should come from some deeper theory. Instead, we shall use present constraints to restrict the range of values 
for M G ,p,2- 

Having seen how one can obtain very small m v , the next question would be: How small can one allow m v to be 
if one takes into account the neutrino oscillation data? First of all, atmospheric neutrino oscillation data gives a 
difference of mass squared Am 2 w 10~ 3 eV 2 while solar neutrino oscillation data gave Am 2 w 10~ 5 eV 2 (MSW) or 
lCP 10 ^ 2 (vacuum). In anticipation of new data, the LSND results are not taken into account in our rough estimation 
of various mass scales. Without any need for a specific model, one can say that the atmospheric data implies that at 
least one of the three neutrinos should have a mass of at least 3 x 10~ 2 eV, while the solar data implies that at least one 
of the remaining two should have a mass of at least 3 x 10 _3 eV (MSW) or 10~ 5 eV (vacuum). As we have seen above, 
the 4th neutrino can be quite heavy. For the sake of argument, let us assume here that its mass is approximately 
100 GeV. Since our three light neutrinos are practically degenerate -a lifting of which will be discussed below, the 
atmospheric data alone constrains R to be greater than approximately 10~ 14 . This in turn constrains M 2 < 10 12 (in 
units of Mp) for the case M 2 > Mq, or Mq < 10 12 for the reverse case. Notice that this rough estimate is only for 
illustration purpose. 

There is however one interesting piece of information which could be quite interesting, phenomenologically speaking: 
the presence of vector-like fermions which carry weak quatum numbers and which could be relatively "light" . These 
are the fermions F with mass Mp as indicated above. Let us recall from the above discussions that Mg,p,2 are all 
expressed in units of Mp which itself could take on any value, even a few hundreds of GeV. The sole restriction will be 
from experimental constraints, a subject to which we shall come back below. Furthermore, we can see from the results 
that the mass of the pseudo-NG bosons can also be "low" as well (Fig.l) which could provide a further experimental 
clue. 

We now turn to an important issue: the lifting of the mass degeneracy of the light neutrinos. The analysis presented 
below will reveal quite interesting implications such as the correlation between the actual values of the masses and 
Am 2 , which can have a profound cosmological consequence. For neutrino masses which are large enough to provide 
part of HDM, the MSW solution of the solar neutrino problem is preferred. If the vacuum solution turns out to be 
the correct one, the neutrino masses will be much too light in our scenario to play a role in HDM. 

We shall divide the discussion presented below into two parts. First we analyze the case when there is no mixing 
between the 4th neutrino and the lighter three. It will be seen that an interesting feature emerges: Am^ ss Am 21 -a 
quasi-symmetric splitting. (Am^ is of the same order.) This phenomenon could be called a mass splitting quasi- 
degeneracy. Of course, solar and atmospheric neutrino data suggest otherwise. Next, we will show how this mass 
splitting quasi-degeneracy can be lifted, suggesting- at least in our scenario-the presence of a 4th neutrino. 

In what follows, we will neglect any possible CP phase in the neutrino mass matrix since we will be concerned only 
with Am 2 and present data on neutrino oscillations are not sensitive to the presence of such a phase. In addition, 
we shall concentrate in the next two subsections only on Am 2 . A full comparison with the data will necessitate the 
inclusion of the leptonic "CKM" angles coming from Vl = U\U V - In the two subsections presented below, we shall see 
what U u might look like. To complete the discussion, we shall use a model for Ui in order to make some statements 
about the size of the mixing angles. The subject of the charged lepton mass matrix itself will be dealt with in a 
subsequent publication. 

B. Neutrino mass matrix I: What if there is no mixing between the 4th and the lighter three neutrinos? 

The 4x4 neutrino mass matrix obtained at this point is purely diagonal. We would like to examine how mass 
mixing might arise. In particular, we would like to lift the degeneracy of the three light neutrinos. In this section we 
will concentrate on the scenario where there is mass mixing only among the three light neutrinos. We will show that, 
in this scenario, Am 23 ss Am^. If this were experimentally the case, it would be hard to detect the influence of the 
4th neutrino since it does not mix with the other three. Since the atmostpheric and solar data appear to point to 
Am| 3 S> Am| 1; we will present in the next section what can be done in order to be in agreement with the data. It 
turns out that this can be accomplished if one introduces a mixing with the 4th neutrino. This implies that, at least 
in our model, Am 23 3> Am 21 implies the existence of a 4th neutrino, and hence a 4th generation. 
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The degeneracy of the three light neutrinos at this level comes from the fact that there is a remaining global 5*0(3) 
symmetry which manifests itself through the equality of the masses of the family gauge bosons (Ma) as well as those 
of the pseudo-NG bosons (Mp). It is then clear that one needs to break that remaining global symmetry in order 
to remove the degeneracy of the light neutrino masses. We would want to do this in such a way as to preserve the 
quasi-degeneracy of the light neutrinos. There are probably several ways to achieve this, and we will present one of 
them here. 

Since we have seen how the diagonal elements of the neutrino mass matrix for the three light neutrinos are obtained 
at the one loop level, it is natural to envision a scenario in which the mixings themselves are obtained at one loop. 
A look at Figs. 1 reveals that the most "straightforward" way to induce mixings at one loop is for (li and/or Repi 
to have mixed couplings, i.e. to both vn and v^j as well as to both -qm and r\^. This could come from mixings 
among flj with different family indices and/or the mixings among Repi. Before getting into the details of what 
kinds of interactions are needed to break the remaining global 5*0(3) symmetry and hence inducing the mixings, it is 
instructional to assume that such a mixing among the boson masses occurs and to write down the Yukawa couplings 
(|l3|,[l4|) in terms of the new boson mass eigenstates. 

Let us first look at the states f2j. As we have discussed earlier, these are the NG bosons which are absorbed by 
the corresponding family gauge bosons. When these NG bosons get mixed, there will be mass mixings among the 
corresponding family gauge bosons. Let us denote the orthogonal matrix which diagonalizes these family gauge bosons 
by Aq . We shall choose the following representation for An ■ 

■S1S2C3 + C1S3 C1S2C3 + S1S3 
An = \ -C2S3 C1C3 + S1S2S3 -C1S2S3 + sic 2 I (24) 
-S1C2 C1C2 

where c and s represent the cosine and sine. If we denote by the logitudinal components of the gauge boson mass 
eigenstates, its relationship with Qi in the unmixed case is given by 

(25) 

y r^3 j \ j 

where Aq is given by 

(C2C3 -C2S3 Si 

-S1S2C3 + C1S3 c 1 c 3 + s 1 s 2 s 3 -sic 2 I (26) 

C1S2C3 + S1S3 -C1S2S3 + S1C2 cic 2 

The masses of the corresponding gauge bosons are now denoted by 

M 2 Gi = M% + 6 1 ; M 2 2 = M% + S 2 ; M 2 3 = M% , (27) 

where 5i t 2 can be positive or negative. Notice that Si^ 2 and the mixing angles shown above are related, i.e. they are 
all derived from the same boson mass matrix. We will show an example of such fact below. 

We can now replace the unprimed states in Eqs.( ]T^Jl4] ) by the primed states using Eq. (p5|). We can then compute 
the one-loop contributions to the elements of the neutrino mass matrix M. v . Let us first look at the contributions to 
the light neutrino masses and mixings coming from the Qi states. The two terms which are crucial for this computation 
are 

GJl cosPQiFa = GJl cos A^^Fr (28) 

and 

G 3 M 2L sin (3^1^ = G3M2L sin Pfy^d? (29) 
In the loop integrations, one will encounter the following propagators: 

fc 2 - M G3 = k 2 — M G ' (30a) 

W = 1,2 _ AJ2 + (1,2 _ J1//-2 Vk2 _ A/T2 V ( 30b ) 















fi 3 / 




V^3 



k 2 -M G1 k 2 -M 2 (k 2 — M G2 )(k 2 — Mg) 
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P - M% 2 k* - Ml + (fc 2 _ Af 2 2)(fc2 _ M 2 ) > (30c) 

With the above remarks in mind, let us proceed to calculate the contributions of f2' to the neutrino mass matrix. 
We shall concentrate first on the 3x3 submatrix of the light neutrino sector. As a prelude to the computation of the 
full submatrix, let us show how two elements are calculated: M.]} and Ail 2 . In these computations, we shall use, as 
an example, the explicit form for An shown in Eq. (24). For the complete calculations of the matrix elements, we 
shall use the notations A^ for An- 

a) In the calculation of the contribution of f2' to M.]}, one combines Eq. (|2^) with Eq. ( p6| ) to get the following 
combination of Q': 

(c 2 c 3 n' 1 -c 2 s 3 n' 2 -s 2 n' 3 )\ (31) 

which gives the following combination of propagators: 

clcKn'^) + cls 2 3 (n' 2 n' 2 ) + s 2 2 (n' 3 n' 3 ) (32) 

Upon using the propagators listed in Eqs.(^o|) in the one-loop integral (Fig.l), one obtains the following replacement 
(the reader is referred to Eq. (um for a comparison): 



32 7T 2 32 7T 2 



■ J ( fi ) ^r^( / ( r! ) + <4<45ii(M G , m G i) + 44^i(M G , m G2 )) (33) 

where (i — 1,2) 



! M F [M 2 {M 2 Mm) ~ M Gi HW)) + M h M h HW)] 
SiI(Ma, Ma) = j^—^i (^-M^M-^MI) °~ ~ " M ^ ^ 

One can see that, in the symmetry limit where Si — > {Mqi — * M G ), 5»J(Mg, M<3i) vanishes identically. 

One interesting remark worth mentioning is the following: In (|3^), the first term contains no mixing angles. 
In fact, the coefficient in front of is c 2 c 2 + c 2 s 3 + s\ = 1, which is the result of An being an orthogonal matrix. 

We do not give the explicit form for because, after taking into account the contribution of Repi, one obtains 
the combination — I(Rep) which is already given by Eq. (17). 

When the boson mass differences, represented by Si, are small compared with M G , another useful form which could 
be used is given by (i = 1, 2) 



SiI(M G , Ma) = -xJ(M G , Xi) (35) 



where 



M 2 M F [-M 2 Al+Xr+Mm))+MG(l+^)] 

= w^w 2 { m^w^TxMim^Mm - {Mf " M2)} ' (36) 

and where 



(37) 



so that 

M G3 = M G ; M G1 = M G (1 + x x ); M 2 G2 = M G (1 + x 2 ). (38) 

Here one could explicitely see the vanishing of SiI(M G , M G i) in the symmetry limit because of the explicit appearance 
of Si on the right-hand side of the equation. 

The other diagonal elements of the neutrino mass matrix can be analogously calculated. One just needs to replace 
the combination of angles in ([32j) with the appropriate ones. 

b) For the 1-2 element, the appropriate combination of propagators is given by 

c 2 c 3 (-sis 2 c 3 + cis 3 )(fi' 1 fi' 1 ) - c 2 s 3 (cic 3 + sis 2 s 3 )(n' 2 n' 2 ) + sis 2 c 2 (tt' 3 D, 3 ) (39) 
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It is now straigthforward to compute M. u ■ It is given by 

Ml 2 (n) = ^^^(c2C3(-S 1 S2C3 + C 1 S3)5 1 I(M G ,Mgi)-C2S3(ciC3 + S 1 S2S3)S 2 I(Mg,Mg2)), (40) 
62 7T Z 

where we have the appearance of the same ^/(Afg, Ma). Notice that A4l 2 ((l) denotes the contribution coming from 
f2 only. The full clement will also include the contribution coming from Rep. 

Notice that the term I (ft) is not present in (f40|). Again this is due to the orthogonality of ^4^. The coefficient 
appearing in front of I {ft) is C2C3(— S1S2C3 + C1S3) — €283(0103 + 5152*3) + S\s 2 c 2 = 0. The orthogonality of Aq implies 
that the product of any two columns is equal to zero. As a result we can see that, in the symmetry limit, M X2 
vanishes identically. This applies to all the other off-diagonal elements. 

In order to complete the computation of the matrix elements (including the 1-1 and 1-2 elements), one has to say 
something about the contributions coming from the pseudo-NG bosons themselves. One might imagine that the same 
mechanism which breaks the global SO (3) symmetry also induces mixing among the degenerate pseudo-NG bosons. 
We will assume that the same matrix A^ diagonalizes the pseudo-NG boson sector so that, instead of the combination 



of 0^ and Repi used in Eqs. (12) for the NG and pseudo-NG bosons, we shall use A^Q, and A^Rep, where f2 and Rep 
are now column vectors. With these definitions, one simply gets W Rep = W A^ 1 A^Rep . This simple assumption is 
used for two purposes: 1) To reduce the number of arbitrary parameters; 2) To see how far one can go with it before 
one needs to modify it. With this assumption, the mass splitting among the pseudo-NG bosons are given as in Eq. 
(|38|), namely 

M 2 3 = M 2 ; Mp 1 = Mp(l + Xl );M 2 P2 = M 2 {1 + x 2 ), (41) 

with the same Xi as for the gauge boson masses. Furthermore, the mixing angles are the same as above. The 
contributions of the "rotated" pseudo-NG bosons to the neutrino mass matrix elements will therefore be accompanied 
by a factor — s 'i^ 2 ^ , just as in Eq. (|3^ 



32tt 2 

As mentioned above, in the full computation of the matrix elements, we shall use, for convenience, Aij to denote 
the matrix elements of Aq instead of the representation of Eq. (|24|). One should then recall that, because Aq is an 
orthogonal matrix, one has: Y]j A 2 ^ = 1 and /~2 k A^A^j — 0. The form of the neutrino mass matrix elements will 
make use of these properties, just as we have done above. 

With the above remarks in mind, the full 4x4 neutrino mass matrix is now given by: 



M v /m N 



/ mn "M2 "ii3 

mi2 to 2 2 m 2 3 

m i3 TO23 ^33 

V 1 



(42) 



where mjy is the mass of the 4th generation neutrino shown in Eq. ([?]). In Eq. (f42|), we have ignored the tiny one- loop 
contribution to 77744 = 1, in particular when there is no mixing between the 4th neutrino and the lighter three. As 
we shall see later on, it can be ignored even if there is mixing, the reason being the fact that my, i,j = 1,2, 3, are so 
much smaller than 77744 = 1. A change of 77744 to a value slightly less than or greater than one will not significantly 
affect the eigenvalues, as we shall see in the numerical examples below. 
With 

AI(G,P,Xi) = I{M G ,Xi) - I(M P , Xi ), (43) 
where I(Mp,Xi) is given by Eq. ( fj6| ) with the substitution Mq —> Mp, one obtains for my: 

m n = ^^{A/(G, P) - A 2 nXl AI(G, P, Xl )- A 2 12 x 2 M(G, P, x 2 )} (44a) 

32 7T Z 

77122 = 1U^~ {AI{G > P) ~ A 2i*iMG, p,xi)- A 2 22 x 2 AI(G, P, x 2 )} (44b) 



™33 = ^^{A/(G, P) - A 2 3lXl AI(G P,xi) — A 2 32 x 2 AI(G, P, x 2 )} (44c) 

61 7T 



mia - -^^-{A 11 A 2lXl AI(G, P >Xl ) + A 12 A 22 x 2 AI(G, P, x 2 )} (44d) 
62 ir z 
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mis = - S -^^-{A u A 3lXl AI{G, P, Xl ) + A n A 32 x 2 AI(G, P, x 2 )} (44e) 

52 IT* 



m 23 = - Ti 2 ? {A 2 iA 3lXl AI{G, P, xi) + A 22 A 32 x 2 AI{G, P, x 2 )} (44f) 

52 IT 



where Ay denote the matrix elements of Aq, as mentioned above, and where AI(G,P) was already defined in Eq. 

A few remarks are in order here. First, one can see that, in the limit Xi — > 0, M. u reduces to a diagonal matrix with 
three equal diagonal elements: h ' 3 ^f AI(G, P). Secondly, apart from various mixing angles, the off-diagonal elements 
depend on results of loop integrals, AI(G, P, x^) which, in turns, depend on the same parameters as the ones that enter 
the loop integrals of the diagonal elements in the unbroken case, AI(G, P). The ratio Rj = AI(G, P, Xi)/AI(G, P) is 
plotted in Figs. (6,7), for two values of the parameter x, as a function of M 2 in the similar manner to Fig. 2-5. (The 
two values of x were chosen for the purpose of illustration and to coincide with the two examples given below.) It 
can be seen that the ratio Ri is at most of 0(1CP2), even for x as large as 0.5. Therefore, in our model, a small mass 
splitting in the scalar and gauge sectors results in a scenario with almost degenerate light neutrinos. The difference 
of the mass squared, Am 2 , depends, however, on the size of the off-diagonal elements. To see how it actually works, 
a simple model of mixings will be presented below along with some numerical examples. 

We starts out with a very simplistic model of mixing and try to see how far one can go. It is: 

A4q P — Mq p I b 1 I (45) 

V° i/ 

where & is a small parameter less than unity. This simple model has the merit of elucidating the points that we 
have made above. (An extension of this model, showing similar results, will be discussed below.) The above mass 
mixing ( |45"1) could come, for example, from a term in the Lagrangian of the form: \${(p, a p' a )(Q,@ p'^) + (p a p' a )(p l3 p'p))- 
Assuming < p' >= (v 1 , 0, 0, 0), < p" >— (0, v", 0, 0), with v'>" <§; V, V , one can obtain the above mass mixing matrix. 




It is easy to see that the eigenvalues of ( |45[ ) are 



Mgx,pi = M 2 . P (1 + b),M% 2>P2 = Mq p(1 - b),M% s>PS = Mq P , (46) 
Aq as discussed above is now given by 




^ = -4* 4, (47) 



Now we can make the following identifications: x\ = b, x 2 — —b. The various angles are given in Aq. The matrix 
elements of the neutrino mass matrix are now fairly simple: 

ron = 8 -^P^{AI(G, P) - h(AI(G, P, b) - AI(G, P, -&))} (48a) 

61 IT I 



m 22 = ^M{A/(G, P) - ^b(AI(G, P, b) - AI(G, P, -&))} (48b) 
62 tt z 2 



m 33 = S ^-{AI(G,P)} (48c) 

62 7T Z 



mia = - S ^^-{h(AI(G, P, b) + AI(G, P, -b))} (48d) 
62 it 2 

mi 3 = 0, (48e) 

m 23 = 0. (48f) 
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The above matrix elements are surprisingly easy to handle. When they are substituted into Eq. (|42|). one obtains 
straightforwardly the following mass eigenvalues: 



mi = mjv 



sin(2/3) 

32 7T 2 



{AI(G,P)-bAI(G,P,b)}, 



m 2 = m N - 



sin(2/3) 

32 7T 2 



A/(G,P), 



m 3 = m N ^^{AI(G, P) + bAI(G,P, -b)}, 



= rriN 

The matrix which diagonalizes the above neutrino mass matrix is simply 



Uy 



f 75 7= °\ 

10 

, -A 71 . 

V o ooi/ 



One obtains the following mass splittings: 



2 2 / 

m 3 — m 2 = (mjv 



!^M) 2 (26A/(G, P)A/(G, P, -6) + (bAI(G, P, -b)) 2 ), 

32 7T Z 



(49a) 

(49b) 

(49c) 
(49d) 



(50) 



(51a) 



m\ -mj = (m N -^) 2 (2bAI(G,P)AI(G,P,b) + (bAI(G, P, b)f). (51b) 

32 7T 

In general, AI(G,P, Xi) < AI(G,P), and combined with the fact that 6 < 1, one has (bAI(G,P,-borb)) 2 <C 
2bAI(G,P)AI(G,P,-borb). One can then neglect the last terms in Eq.@. Numerically, one has AI(G,P,b) m 
A/(G, P, —6). This implies that m 2 — m\ w m| — m 2 , a quasi-degenerate mass splitting. This holds for any value of 
b. Solar and atmospheric data suggest otherwise. This necessitates the lifting of this quasi-degeneracy of the mass 
splitting. To do this, we need to invoke some kind of mixing between the 4th neutrino and the lighter three. In an 
indirect way, the disparity between Am| 3 and Am^ indicates- in our model- the influence of a 4th generation. Before 
discussing this issue which will be presented in the next section, let us illustrate numerically a few examples of the 
quasi-degenerate case. 

First, a few useful points are in order here. Since m 2 = m n 8 I^^P AI(G, P) , one can rewrite the above equations 
( |5l| ) as (neglecting the last terms on the right-hand side) 

ml ~m\= m 2 {m N 2b S1 ^ 2 ^ AI(G, P, -b)), (52a) 



m 2 2 -m 2 = m 2 (m N 2b-y^AI(G,P,b)). (52b) 
32 n z 

For a fixed value of m 2 , the size of the mass splitting, Am 2 , depends on the size of the factor 
mjv(26) (sin(2/3)/32 7r 2 )A/(G, P, — bor b). At first glance, it appears that one can obtain Am 2 to be as small as 
one wants with the appropriate choice of b. Although it is true that it can be so, we will show that, Am 2 can also 
be very small (< l0~ lo eV 2 ), even when b w 1. This depends on how large the masses of some of the particles 
participating in the loop diagrams are. As a result, by limiting Am 2 > 10 _10 ey 2 , one puts a constraint on those 
masses. 

In Fig. 8, we present the "median" mass m 2 as a function of M 2 and Mq for a given Mp (as presented in Figs. 
(2-5)). The mass is given in units of (m.Ar/100 GeV). Similarly, we present in Figs. (9, 10) m§ — m\ and m\ — m\ as 
a function of the same masses, but also for a given value of b. The results are expressed in units of (m^ /100 GeV) 2 . 
For a more streamlined presentation of the results, we shall limit ourselves to the case m 2 < 1.67 eV, coming from the 
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suggestion that the sum of neutrino masses lies between 4 and 5 eV in order to form a component of HDM. Similarly, 
we shall restrict Am 2 < 1 eV 2 . In our model, for a given value of b, and Am 2 are correlated as one can see from 
Fig. (8,9,10). 

Three major remarks are in order here. 1) One can see from Figs. (9, 10) the quasi-degeneracy of the mass splitting 
in this particular scenario. (In the next section, we shall see how one can lift that degeneracy.) 2) One can also see 
from Figs. (8, 9, 10) that, were the vacuum solution to the solar neutrino problem favored, i.e. Am 2 sa 10~ 10 ey 2 , 
the median value ui2 will always be less than 0.1 eV. (The lifting of the mass splitting degeneracy to satisfy the 
atmospheric neutrino data will not change this conclusion.) This simply means that, at least in this model, the solar 
vacuum solution is incompatible with the light neutrinos being significant components of HDM. 3) Also from Figs. 
(8, 9, 10), it can be seen that the MSW solution, Am 2 w 10~ 5 eV 2 , can correspond to values of mi larger than 1 
eV. (Again, the lifting of the mass splitting degeneracy to satisfy the atmospheric neutrino data will not change this 
conclusion.) So, in our scenario, the MSW solution is compatible with the light neutrinos being significant components 
of HDM while the vacuum solution is not. This is a very specific prediction of this model. 

The above discussion leaves out the question of the size of the mixing angles. As mentioned above, we have already 
fixed the neutrino mixing matrix U v , as given by Eq. (|50|). To complete the task, one has to model the charged lepton 
mixing matrix [//. This is something that we shall do in the last section. We wish however to reemphasize the main 
result of this section: the values of Am 2 are independent of Ui. As one can see from Fig. (9, 10), Am 2 depends only 
on the various masses and on the parameter b, regardless of Ui. As a consequence, the large angle or small angle 
solutions as deduced from the data basically constrains, in our scenario, the matrix TJ\ (JJ V being already fixed). 

To finish the discussion of this section, we wish to present another form for the boson mass matrix, namely 



M 



G,P 




The mass eigenvalues are 

Ml hPl = Mq P (1 + V2b),M 2 2 ^ P2 = M 2 G >P (1 - \2b:.\l 2 :u , :i = M 2 j 

An is now given by 



An = 




(53) 



(54) 



(55) 



V y/S 

It is now straightforward to see that the neutrino mass matrix elements are 



mn = m 33 = 



sin(2/3) 



{AJ(G, P) - — = (AI(G, P, b) - AI(G, P, -&))} 



32tt 2 l y 7 ' 2V2 



(56) 



m22 = SJ ^ L{AI{G ' p) ~ 7! (A/(G ' F ' b) ~ A/(G ' p ' ~ b))} (57) 

mia = m 23 = -^^{h(AI(G,P,b) + AI(G,P,-b))} (58) 

mia = -^^{^=(AJ(G, P, b) - AI(G, P, -&))}, (59) 

The eigenvalues are now simply given by 

mi = m Ar ^M{A/(G, P) - V2bAI(G, P, b)}, (60a) 
32 71" 2 

m 2 =m N ^^AI(G,P) 1 (60b) 
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m 3 = m N ^^l{AI(G, P) + V26A/(G, P, -b)}, (60c) 

32 7T Z 

7774 = itin (60d) 

These masses have exactly the same form as those of Eq.(|49"|), except for the factor of ^/2b instead of b. The matrix 
U u which diagonalizes the above matrix is exactly the same as in Eq. (|50| ) . Furthermore, m\ — mf, and m 2 — m\ are 
of the same form as Eqs. (f5l|), with the following replacement in Eqs. (plf): b — > 6' = The analysis which follows 
is exactly the same as the one presented above. 

One can envision various scenarios for the boson mass matrices, but it is certainly beyond the scope of this paper. 
To make things more complicated than the simple assumption (h5J) does not appear to add much to the discussion. 
Although it might be possible that a more involved ansatz than fl4q) could lead to the lifting of the mass splitting 
"quasi-degeneracy" , we have not succeeded in finding it. For this reason, we now turn our attention to the more 
appealing scenario, at least within our model: the mixing between the 4th neutrino and the rest. 



C. Neutrino mass matrix II: Mixing between the 4th and the lighter three neutrinos 

We have seen above that the simple ansatz for the boson mass matrices (^) leads to a situation in which the mass 
splittings are quasi-degenerate. This, of course, is in contradiction with the data. In this model, in order to lift that 
quasi-degeneracy, one needs a mixing between the 4th neutrino and at least one of the lighter three. To get a feel for 
what might be needed, we shall first present a few numerical examples. Based on these examples, we shall attempt 
to give a theoretical basis for these numerical examples. 

As an example, we shall choose a specific value for the parameter b and for the masses M 2 , Mq, Mp and Mp which 
enter the loop integrals for the neutrino masses. This will fix a definite value for the matrix elements of the neutrino 
mass matrix. As we have already discussed earlier, the integrals depend only on the ratio of the above masses. We 
will present two examples for the purpose of comparison. We shall see the reasons why we wish to do so below. 

1) First Example: 

We shall set (in units of Mp): Mp = 1, M P = 5, M G = 10 6 , M 2 = 2.5 x 10 9 . For b, we shall choose: b = 0.035. 
(A smaller value of b will give a smaller mass splitting.) The reason for this choice (other choices are equally valid) is 
the fact that it will give a typical mass of approximately 1.5 eV and a desired mass splitting. Putting these values 
into the expressions for the integrals as given by Eq. ([56]) , we obtain the following neutrino mass matrix, where mjv is 
assumed to be 100 GeV for convenience: 



M v = (lOOGeV) 



/ -1.579332216 x lO" 11 .8697647852 x 10 -1 ' 

.8697647852 x 10" 17 -1.579332216 x 10~ n 

-1.579332184x 10" 

V 



(61) 



Notice that the above matrix has no mixing between the 4th neutrino and the lighter three. The eigenvalues are just: 
|mi| = 1.579331346eV; \m 2 \ = 1.579332184eV; |m 3 | = 1.579333086eV; |m 4 | = lOOGe^. (62) 
As we have discussed in the previous section, this gives a quasi-degenerate mass splitting, namely 

Am 2 32 = 1.601195367 x lO^eF 2 , (63) 



Am|i = 1.535757079 x lO^e^ 2 , (64) 

where Am^ = m 2 — mf. 

Let us now assume that the mixing with the 4th neutrino is non-zero. We start out with the simplest assumption, 
namely one in which only the 3rd neutrino mixes with the 4th one. This means that 777,34 and 77743 are both non- 
vanishing. If we wish to have m\ — m\ ss 10 -3 eV 2 as suggested by the atmospheric neutrino data, it turns out that 
m.34 and TO43 cannot be too small nor too large, being of order 10~ 7 totv. Notice that 77134 and 77143 do not have to be 
equal. We shall see how it might be possible to obtain such a number. Let us first see how it works from a numerical 
viewpoint. 

To guide our understanding of how things work, let us notice that, by adding 777,34 and 77743 to M. v above, one changes 
only one of the three light mass eigenvalues, leaving the other two the same. Now the two unchanged eigenvalues will 
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be the ones that fix one of the two mass splittings, Am 2 . For convenience, we shall choose the Am 2 corresponding 
to the unmodified mass eigenvalues as the one which corresponds to the solar neutrino problem. As we have learned 
from the above analysis in Section (I11B), if one chooses the MSW solution, then one can find masses which are 
large enough for HDM, while, if the vacuum solution is chosen, the masses will be too small to form any significant 
component of HDM. For the numerical example given here, we shall choose the MSW solution as shown above. For 
77134 and TO43, we shall first choose a symmetric case (there is no particular reason for this being so) as an example. 
We have 



(lOOGeV) 



( -1.579332216 x 1CT 11 
.8697647852 x 1CT 17 





.8697647852 x 1(T 17 
-1.579332216 x 1CT 11 










^1.579332184 x 
.8 x 1(T 7 



10 



-11 






x 10" 
1 



(65) 



The eigenvalues are 



\mi\ 



1.579331346eF; \m 2 \ = 1.579333086eV; \m 3 \ = 1.579972184eV; m 4 = lOOGeV. 



(66) 



We then get 



Ami 



2.02 x 10~ 6 eV 2 



(67) 



Am 2 -! = 5.497 x 10~ 6 eV 2 . 
The matrix which diagonalizes the mass matrix is 

( 



1 


1 










f 






















-1 


0.8 x 10 





. 


8 x 10" 


" 7 1 



(68) 



(69) 



/ 



Two remarks are in order here. Firstly, from the values of the light neutrino masses, one obtains yV. , |mj| ~ 4.7eV, 
which is in the range of mass for HDM. Secondly, Eq. ( |67[ ) corresponds to the best fit for the atmospheric neutrino 
data, while Eq. ( |6q ) corresponds to the best fit for the (small angle) MSW solution to the solar neutrino data. One 
word of caution: this is not a prediction because we chose the masses (Mg, M 2 , etc..) in such a way as to "reproduce" 
the experimental results. It nevertheless shows a dynamical basis for these numbers. Also, for nothing more than a 
numerical example, the values of 7773443 were chosen arbitrarily in order to have the desired mass splitting. How to 
justify these values is the subject to be discussed below. 

The next numerical example deals with the case when 77734 7^ 77143. in doing the analysis, we notice that it does not 
matter whether 77134 is greater than 77743 or the other way around. One obtains the same result either way. We shall 
require that Am| 2 (e^ 2 ) = 10~ 3 — 10~ 2 . It turns out that 77734 and 77143 can range (in units of mjv) only between 
approximately 0.4 x 10~ 6 and 0.8 x 10~ 8 . To be explicit, one has 



M v = (lOOGeV) 



/ -1.579332216 x 10~ n 
.8697647852 x 10~ 17 


V 



gives A777 2 2 (el /2 ) ~ 10 2 , while 



M v = (lOOGeV) 



/ -1.579332216 x 10" 11 
.8697647852 x 10~ 17 


V 



.8697647852 x 10~ 17 
-1.579332216 x 10~ n 







.8697647852 x 10~ 17 
-1.579332216 x 10~ n 










4.579332184 x 10" 
.8 x 10~ 7 






-1.579332184 x 10 
.8 x 10~ 8 




(70) 



-11 




(71) 



gives A777§ 2 (ey 2 ) w 10 -3 . Notice that Am|i stays the same. The above numerical results show that 77134 can differ 
from 77743 by a large factor (50 in this case) while keeping A77732 within the desired range. 
2) Second Example: 

In this example, we choose (in units of Mp): Mp = 1, Mp = 5, Mq = 10 4 , M 2 = 1.2 x 10 9 . For b, we shall choose: 
b = 0.000095. For simplicity, we shall assume, as we have already done above, the following values for 7773443, namely 
77734 = 77743 = 0.8 x 10 _7 (100GeV^). The mass matrix is now 
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1.382258467 x 10" 11 .981382953 x 10" 17 \ 

., _ | .981382953 x 10" 17 1.382258467 x KT 11 

' " ~ ( ' I n n 1.382258467 x 10- 11 .8 x 10- 7 

.8 x 10~ 7 1 / 








The eigenvalues are 

mi = 1.382259448eV;m 2 = 1.382257486eV; m 3 = 1.381618467eV;m 4 = lOOGe^, 
with the corresponding diagonalization matrix given by 



(72) 



(73) 



( 72 vS 



The mass splittings are 



\ 




i i_ 

0-1 0.8xl0~ 7 
\ .8 x 10~ 7 



(74) 



1 / 



lAm^l = 1.77 x 10- 3 eV 2 , 



(75) 



5.42 x 10 _t W 2 



(76) 



The above two examples are chosen sololy for illustration. Other values of Am 2 are possible with different choices 
of various masses (Mq, M2, etc..) and/or the parameter b. 

Before turning to the discussion on the possible origins of 77134,43, let us briefly discuss the "tiny" one- loop contribu- 
tion to 77144, namely 61714 as given by Eq. (pcf). One might wonder how it would affect the light mass eigenvalues. It 
turns out however that, as long as <5t77,4 <C 1 (which is the case in this paper), it does not matter what value it takes. It is 
easy to see how. A 2x2 matrix of the form (a, c; c, b), where a, c b, has as eigenvalues: b+c 2 /6+(l/4)a 2 /b)+0(c A , a 4 ) 
and a — c 2 /b — (l/4)a 2 /6) + 0(c 4 , a 4 ). One can see that, for the smaller eigenvalue, a small change in b affects very 
little its value. As an example, we put 0.99 instead of 1 in Eq. (). We obtain Am 2 2 (eV 2 ) « 1.02 x 10~ 3 instead 
of 1.06 x 10~ 3 (for 1). If we put 1.1 instead of 1, we obtain Am| 2 (eV r2 ) w 0.923 x 10~ 3 . Considering the kind of 
accuracy that one has at the present time, this is completely irrelevant. 

There are probably several scenarios for calculating 77734,43. However, considering the fact that the present exper- 
imental status is not accurate enough for a detailed model, we will present below a more or less "generic" scenario 
which will show how one can obtain 77734.43 of the right order of magnitude. 

What might be the origin of 77734,43? It might be obvious up until now that the vacuum expectation values of 
f2 and p shown in Subsection (III. A) cannot generate such a mixing. One needs at least one additional scalar with 
a non-vanishing vacuum expectation value along the 3rd direction. Let that field be f2' and let us assume that 
< D.' >— (0, 0, v, 0). Let us also assume that there are couplings of the type: 



(77) 



where, for convenience, we have omitted the SU(2) u r index in p. With the above couplings, one can construct 
diagrams for 77734 and 77743 as shown in Fig. 11. 

We shall denote the masses of H4 and hi by Mh a and M^ 4 respectively. Let us define the following quantities: 



AM 2 (G, Hi) =M%- M- 



(78a) 



AM 2 (G, hi) = Mq — M, 



h 4 3 



AM 2 (P,Hi) = M 2 P -M 2 Hi , 



AM 2 (P,h 4 ) = M 2 P - Ml 



(78b) 
(78c) 
(78d) 



From Fig. 8, we obtain: 
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X 3i u vM F M 2u 2 2 A/(G,P 4 ) 2 2 AI(G,h 4 ) 2 2 A/(P,P 4 ) 



"W"»" - ^)^ S « AM2(G ^ 4) + ^am^g^) - ^AM 2 (P,P 4 ) 

a 2 AI(G,fo 4 ) v , , 

^ AM 2 (G, hi) i79 ) 

, , A34 v WfM2 v 2 a A/(G,ff 4 ) 2 2 AJ(G,ft 4 ) , 2 2 AJ(P,ff 4 ) 
"W«* = (^)(— ^)(^4 AM2((? ^ 4) + ^^ AM 2 (G; ^ 4) + c ^ AM 2 (P,ff 4 ) 

^ tt AM 2 (G,ft 4 ) j V ; 

where c and s stand for cos and sin, and A/(G, H4) and the other similar quantities in Eq. j79| ) are given by Eq. 
(|l8l), with the substitution of the appropriate masses taken into account. 

As one can see from the above equations, the expressions appear rather complicated at first look. However, one 
can make an estimate as to which term in 17x34 and m 4 3 is the most important. Each term in Eqs. ( |79| ) is of 
the form: X(v/V)(Mp/M2)(M 2 /AM 2 )AI(mixing angles), where A stands for A3 4 , 4 3. First, we have seen from the 
above numerical analysis that, if we wish to have a mass of 0(1-2 eV), then M F /M 2 « 1(T 9 . It is reasonable to 
assume that X(v/V) x (mixing angles) < 1. If one of the terms in Eq. ( f79| ) were to be the dominant one and that 
"134,43 ~ 10~ 7 , then one should have (M 2 /AM 2 )AJ > 10 2 . Let us first look at the (G;H,h) contribution. Assuming 
that M HiM < Mq so that (Mf/AM 2 ) w M^/M G , it turns out numerically that (Mf/M^)A7 is always less than 
~ 10. For M Hi ^ hi > M G , AI is larger in value than the previous case, but then with (M|/AM 2 ) sa Mf/Mjy h , one 
will again have (M|/M^- hi )AI less than 10 2 . Taking into account the actual calculation of 77134,43 which includes 
mixing angles and various factors, the (G; H, h) would be too small to actually affect the mass splittings. This leaves 
us with the contribution coming from (P; H, h). Here, as we have done above, we will set Mp = 5 in units of Mp. 
There are several possibilities that one can explore. We will present here one of such possibilities. The main purpose 
will be to show that, under reasonable assumptions, one can obtain the desired order of magnitude for 77i3 4j4 3. In 
addition, one would like to see phenomenological implications coming from such a scenario- something extra other 
than just a mass matrix. 

Let us assume that, by an appropriate choice of parameters in the Higgs potential, one has M# 4 to be of 0(Mq), 
and that M^ 4 -C M 2 . Furthermore, let us assume that one also has j3 rs a. Although it is not really necessary, let 
us further assume that A3 4 ~ A 4 3. Now numerically, (Mf/AM 2 ) A/ < 10 2 when one of the masses in AM 2 is much 
larger than the other one and not too much different from M 2 . This is just the case for Mjj 4 = 0(Mq) 3> Mp. Under 
these assumptions, we are left with the (P; h) contribution. In this case, one has 11134 w 77743. So we get 

|m 34 | « |m 43 | « ^ A 34^I M 2 M| M2 |AM 2 (G,fa 4 ) S 2 c 2 . (80) 

Typically, AM 2 (G,h 4 ) = O(10~ 7 - 10~ n ). In most of our examples, M F /M 2 ~ 10~ 9 . So one would expect 
(M F /M 2 )AM 2 (G : h 4 ) ~ 10 -16 - 10~ 20 . If we wish m 34 w m 43 w m N .8 x 10~ 7 , for example, the other factors 

have to be sufficiently large. First, the ratio | M i^ 2 M i | can be rather large if Mh 4 is small compared with M 2 . 

p /14 

Secondly, even if the previous ratio can be large, it can still be offset by s^c 2 . Let us recall from Eq. ( [Hf ) that 
tan/3 g 2 (M F /M 2 )(M^/M G ) « 2 1O -9 (M|/M£). Therefore the angle can be very small if Mq is too "close" in 
mass to M 2 . A numerical investigation reveals that, if one wants to have a mass of 0(leV) and, at the same time, a 
large enough angle, Mq can be relatively "low" (~ 10 4 in units of Mp). (This would imply that the scale of family 
symmetry could be a few thousands of TeV if Mp is a few hundred GeV's.) We now give a couple of numerical 
estimates. We shall take the Second Example as a prototype. There one can calculate the factor s^c 2 to be « 0.134. 
1) For M/j 4 = 100 with all other masses being the same as those of the Second Example, we obtain 

to 34 ss m 43 w toatA3 4 -^: x 4.7 x 10~ 7 . (81) 

If we wish m3 4 » m 4 3 « mjy.8 x 10~ 7 , then A3 4 -^ « 0.17. So one could either have A3 4 s» .2 and v s» V, or some 
other combination. 2)For M/ l4 = 10, we have 

m 34 Rim 4 3Rim J vA34^xl.4xl0" 5 , (82) 
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which would imply ~ 0.006- a reasonable constraint. 

It turns out that the cases with M^ 4 > 1000 (in units of Mp) do not work because then the mass ratios are not large 
enough to compensate for the smallness of the integrals. It is interesting that one can have scenarios where hi is light 
enough (i.e. not too much heavier than F)- a feature which could have interesting phcnomcnological implications. 



D. Oscillation Angles 

To discuss the neutrino oscillation angles, one needs to give the leptonic "CKM" matrix, namely Vl — U\\J V . It 
is beyond the scope of this paper to discuss the charged lepton sector, and hence Ui. This will be the subject of the 
following publication. However, we can give an example of Ui by adopting, at least for this paper, a simple model 



of charged lepton masses of Ref. ( 15 1), which is a phcnomcnological model based on a generalization to the lepton 
sector of the "democratic mass" ansatz of the quark sector. The reason why we use, as an example, Ref. ( |lf|) is 
because the matrix which diagonalizes the neutrino mass matrix, U v , is identical to the 3x3 submatrix of our Eq. 
( p0| ) (apart from a difference in in the overall sign) , namely 

U^- ± -4, (83) 




Although Ref. ( fllql ) discussed an ansatz for three generations, we will use it here because the mixing with the 4th 
generation is not relevant for the oscillation angles we are interested in. (It was relevant for the mass splitting.) So, 
basically, we will be using only the phenomenological ansatz for the charged lepton mass matrix of Ref. ( p5[). In 
fact, we will only use the matrix which diagonalizes that mass matrix. 
The 3x3 leptonic "CKM" matrix written down by Ref. ( §]|) is 

/ _J -(l/\/3V"*e/nv (2/V6) y /m e /m ll \ 

Vi = {ABrfU v « y/me/mn l/y/3 -2/y/E (84) 

V 2/V6 1/V3 / 

where ABi is the matrix which diagonalizes the charged lepton mass matrix, U v is given above, and m e and m M are 
the electron and muon masses respectively. Now, the probability for u e — ► v fl is 

P(ve -» vy) « 2(V? 1 Vi 1 + V? 2 Vl - sin 2 (1.27Am 2 2 L/£:), (85) 

where the usual notation sin 2 (20i2) is simply the coefficient of sin 2 (1 .27 Ato 2 2 L / E) . Similarly 

P(v„ -» Vt) » 4y 2 2 3 y 3 2 sin 2 (1.27Am 2 3 L/i?), (86) 

with sin 2 (26*23) being the coefficient of sin 2 (1.27Am 2 3 L/_E). Putting in the values of m e and m M to evaluate the 
matrix elements of Vi , one readily obtains 

sin 2 (26>i 2 ) w 6.5 x 10~ 3 ; sin 2 (26i 23 ) w 0.89. (87) 

These results correspond to the small angle MSW solution, and to the large angle atmospheric solution respectively. 
This is consistent with the best fit for the two neutrino oscillation problems. 

The above results should be viewed with caution. The small angle MSW solution given above, as well as the large 
angle solution for the atmospheric oscillation, depends on the charged lepton sector - the neutrino sector diagonalizatin 
matrix being already fixed by Eq. d5(i|). One can easily imagine how these angles can drastically change if the charged 
lepton mass matrix has a different texture. This will be the subject of a subsequent paper where we will examine the 
charged lepton mass matrix in the context of the present model- the basic interaction Lagrangian being already given 
by Eq. (|). 



IV. EPILOGUE 



The above discussions focused entirely on the atmospheric and solar neutrino data. We have left out the LSND |l6| ] 
result for two reasons. Firstly, it is because it might be prudent to wait for future experiments, either to confirm or to 
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refute these results. Secondly, it is because it is extremely hard to incorporate all three experiments simultaneously in 
a "natural" model. In general, one needs to invoke some kind of sterile neutrino that mixes with the lightest neutrino 
to explain the solar data. If this sterile neutrino were to arise from some kind of model, it is rather hard to invent, in 
a "natural" way, a scenario to explain why this sterile neutrino is so light and close in mass to one of the three active 
light neutrinos. 

Let us suppose that the LSND result are verified by future experiments. What does the model presented in this 
paper have to say about a sterile neutrino? Let us remember that rj R — (v R , u R ) is an electroweak singlet. Furthermore 
we have seen that it is which mixes with to give masses to the neutrinos. Its SU(2) UR partner, v R ), remains 
massless, at least within the framework of the preceding sections. Could these be the so-called sterile neutrinos? If 
so, how would they get a mass? How would they mix with the light neutrinos? These are the questions which are 
under investigation. 

We have concentrated in this manuscript on the even option. One might wonder about the odd option and its 
implication on neutrino masses. It is beyond the scope of this paper to investigate this issue, however a preliminary 
investigation of the odd option, with three families and one family singlet rf, appears to indicate that the preferred 
solution for the neutrino masses is that in which there is a hierarchy m\ -C mz -C TO3. 

There are numerous phenomenological consequences to be worked out in subsequent publications. One can, however, 
make one rather solid prediction: neutrinos, being of the Dirac nature, will not give rise to the phenomenon of 
ncutrinolcss double beta decay. Another interesting consequence is the possible existence of "light" (i.e. 200 GeV or 
so) vector-like fermions: F, as well as TeV-scale pseudo NG bosons which carry family and SU(2) UR quatum numbers. 
This will be dealt with in a separate paper. 

Several other phenomenological issues remain to be investigated. For instance, what are the consequences of a 
broken SU(2) uR and what might the cosmological implications of u R s and r]' R be? When SU{2) vR is broken by pf, 
the gauge bosons are expected to acquire a mass of O(V') and can be quite heavy. Since only right-handed neutral 
leptons participate in SU(2) uR interactions, a place where the effects of those gauge bosons might show up is in 
the decays of neutrinos. Without going into detail, it is easy to see that the decay of the light (near-degenerate) 
neutrinos into each other is completely negligible for lack of phase space and for the fact that neutrino masses are 
tiny compared with V' (even if the latter is in the TeV region). This leaves us with the decay of the (heavy) fourth- 
generation neutral lepton TV for which we have JV — > TV + z/j + vi (1) via the exchange of SU(2)„ R gauge bosons, and 
N — > Zr + lj + Vj (2) if mjy < mw or N — ► /J/" + W (3) if m^ > raw- In addition, one could have N — ► E + l^ + Vj 
when mjv > Trip, via the exchange of W. Whether or not mjv is larger or smaller than roj, the relevant decays to 
compare with each other are (1) and (3). To make an estimate, let us assume the the family gauge coupling is about 
the same size as the electroweak coupling (g ~ 0.7). The ratio of the decay widths for (1) and (3) is approximately 
r(l)/r(3) - 7.5 x 10~ A (M w /M d ) 2 {m N /M d ) 2 (l - (M w /m N ) 4 )- 2 X - 2 , where M G represent the mass of the SU{2) VR 
gauge bosons and x represents the mixing cofficient between the 4th neutrino and a light charged lepton. Now let us 
remember that the computation of the neutrino masses does not involve Mq and as a result there appears to be no 
constraint there. However, Mg ~ gV and M G ~ gV, and as a result tan/3 = V'/V ~ Mq /Mq ~ g 2 10~ 9 (M% /M G ). 
In the second example discussed in the previous example, V 1 <~ V (with Mq = 10 4 M^ ) which implies M G ~ Mq- 
Now r(l)/r(3) can also be appreciable if win is close to mw- For example, if Mp ~ 200GeV and mjv ~ 82GeV, 
r(l)/r(3) ~ 1 provided x ~ 10~ 9 . If this were the case, the signal would be quite interesting: a long-lived massive 
neutral lepton whose electroweak decay width is not what it should be. It is certainly beyond the scope of this paper 
to explore numerous phenomenological consequences which might arise from our scenario. 

As for the cosmological consequences of vrs and n R , if they are massless, one should recall our earlier discussion: 
These particles only have family and SU (2) uR gauge interactions (both for v R s and the latter only for n R ). Therefore, 
they cannot influence big-bang nucleosynthesis. One can estimate their decoupling temperatures by comparing the 
interaction rate Ti nt ~ G 2 T 5 , where G 2 ~ 1/(64V^') 4 ), with the Hubble rate H ~ T 2 /m p i. Decoupling occurs when 
Tint < H which gives a temperature of O(10 6 ) GeV if ~ 10 9 GeV for example. After this, their temperature 
would scale like T ~ 1 / R. It is not clear what else they can do except to exist as almost non- interacting relativistic 
relics with an energy density negligible compared with normal matter. At this stage, it is also not clear if they really 
do need to have a mass. The cosmology of these objects is probably worth exploring further. 

Another interesting cosmological subject to explore is the "heaviest" particle in our scenario: The vector-like neutral 
fermion M2 which is singlet under all the listed gauge groups in Eq. (||). M couples to other fermions via||. The decay 
modes obtained from (||) are: M.i R — > (1) and M.2L — ► Pcjf 1 (2). Notice that, in the examples given above for 

the calculations of the neutrino masses, the mass of this fermion is typically M2 ~ 10 9 Afp. So, if Mp ~ 200 GeV (or 
a few hundred GeV), one would then expect the mass of M.i to be around a few times 10 11 GeV. If Mp ~ 1 TeV, 
M.2 would have a mass around 10 12 GeV. The questions that we would like to investigate are: (1) How many A4 e 
are left in the present universe?;(2) Could the decay of the relic M^s manifest itself as ultra high energy cosmic rays 
(UHECR)(with energy exceeding 10 20 eV — 10 11 GeV) whose origins are still unknown? It does appear that the mass 
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of M.2 is in the right energy ballpark. This would be the case of a non-accelerated source of UHECR and is part of 
the "top-down" approach to UHECR 0. For example, M.2R would decay into the longitudinal component of W 
(0 ± ) and FjF. ± would in turn decay into extremely high-energy quarks and leptons. The quarks will hadronize into 
hadrons such as pions which will eventually convert into photons, neutrinos, etc.. 

Last but not least, in the subsequent series of papers, we shall deal with the charged lepton sector and with the 
quark sector. In particular, we shall see how the generalization of Eq. (^|) to the quark sector might yield interesting 
results. 

I would like to thank Vernon Barger, Paul Fishbane and Paul Frampton for reading the manuscript and for useful 
comments. This work is supported in parts by the US Department of Energy under grant No. DE-A505-89ER40518. 



APPENDIX: HIGGS POTENTIAL 



In this Appendix, we shall discuss a simple form of the Higgs potential for the group SO(4) ® SU(2) JJR . For 
simplicity, we shall assume that there is no cross coupling between (f2, p) and the SM Higgs field (f>. (One might 
wonder about the fact that even if the cross coupling were vanishing, it might still be induced through radiative 
corrections. This, however, would be very small in our model.) 

The potential containing f2 and p reads 



v(n, P ) = A x (frn Q - v 2 ) 2 + \ a tf* Pa - v' 2 ) 2 + A 3 [(rro a - v 2 ) - tf a Pa - v' 2 )} 2 
+x 4 m a n a )(p^ Pf3 ) - (srpt,)(siPpf,) 



(Al) 



where < fi >= (0, 0, 0, V) and < p >= (0, 0, 0, V <X> Si), with s± = ( ^ ) . Here, we will assume that f2 is real and p 



is complex. We will be particularly interested in the mass eigenstates resulting from Eq. (Al) 



With n 4 = H 4 + V and p 4 = 



h 4 + V + i<j)± 



The eigenvectors are 



Eq. (Al) gives rise to the following mass matrix for H 4 and h 4 : 

(A2) 



(A! + A 3 )U 2 -\ 3 VV 
-X 3 VV (A 2 + A 3 )U' 2 



H 4 = cos a H 4 + sin a h 4 



(A3a) 



The associated eigenvalues are 



where 



hi = — sin a H 4 + cos a h 4 



m 2 Hi =4(\ 2 + \ 3 )V ,2 m 2 , 



m 2 = 4(A 2 + A 3 )U' 2 m 2 , 



m? 2 



l + a± ^{l-af+Ab 2 



i Ai + A 3 2 
a = {- — — — )tan /3, 

A2 + A 3 



A2 + A 3 



-)tan/3, 



V 

tan/3 = — , 



(A3b) 

(A4a) 
(A4b) 

(A5a) 
(A5b) 
(A5c) 
(A5d) 
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cos a = — ; (A5e) 
y/l + [(l-ml)/b]* 



V' 2 —W 

2A 4 ( _ T/T/ , U ) (A6) 



The mass matrix for 0^ and Repi with i = 1, 2, 3, is 

V' 

-VV V 2 
The eigenvectors are 

Qj = cos/3fij + sin/3i?epi, (A7a) 
Repi = — sm f3 Hi + cos f3 Repi, (A7b) 

The associated eigenvalues are 

TOj 2 = 0, (A8a) 

m fle p = 2A 4 (V 2 + V' 2 ). (A8b) 

Notice that fii are NG Goldstone bosons which are absorbed by some of the S'0(4) gauge bosons. 

Since it is not of immediate relevance to the paper, we will simply quote the masses of the other scalars obtained 



from (Al). Scalars (Pseudo-NG bosons) which have a mass 2X4V 2 : Impi, Rep' i: ltnp\. Goldstone bosons which are 
absorbed by some of the 5*0(4) ® SU{2) VR gauge bosons: Rep' 47 Imp' 4 . Notice that the pseudo-NG boson masses are 
all proportional to A 4 . As a result, their masses tend to zero as A 4 — > 0. 



[1] 
[2] 

[3] 
[4] 

[5] 
[6] 
[7] 
[8] 
[9 
[10 

[11 
[12 
[13 

[i4; 

[15, 
[16 
[17 



e-mail: pqh@virginia.edu . 

Super-Kamiokande collaboration, Phys. Rev. Lett. 81, 1562 (1998). 

See Peter Fisher, Boris Kayser, and Kevin S. McFarland, hep- ph/9906244, for a comprehensive review and a list of 
references. 

M. Gell-Mann, P. Ramond, and R. Slansky in Supergravity, edited by D. Freedman et al. (1979). 

J. C. Pati and A. Salam, Phys. Rev. D 10, 275 (1974); R. N. Mohapatra and J. C. Pati, ibid. 11, 2558 (1975); R. N. 
Mohapatra and G. Sejanovi'c, Phys. Rev. Lett. 44, 912 (1980); Phys. Rev. D 23, 165 (1981), and references therein. 
P. Q. Hung, Phys. Rev. D59, 113008 (1 999). 
See David O. Caldwell, " 



tiep-ph/9910349, for a review of the subject. 



E. Witten, Phys. Lett. B117, 324 (1982). 

J. I. Silva-Marcos, Phys. Rev. D59, 091301 (1999). 



See Jens Erler, hep-ph/990344£ , for a recent review and a list of references therein. 

P. Q. Hung, Phys. Rev. Lett. 80, 3000 (1998). 

P. H. Frampton and P. Q. Hung, Phys. Rev. D 58, 057704 (1998). 



P. H. Frampton, P. Q. Hung, and M. Sher, hep-ph/9903387 , Physics Reports in press. 
Heidelberg-Moscow Collaboration, Phys. Rev. Lett. 83, 41 (1999). 

B. Adeva et al (L3); D. Decamp et al (ALEPH); M. Z. Akraway et al (OPAL); P. Aarinio et al (DELPHI), Phys. Lett. 
B231, 509 (L3), 519 (ALEPH), 530 (OPAL), 539 (DELPHI) (1989). 
M. Fukugita, M. Tanimoto, T. Yanagida, Phys. Rev. D57, 4429 (1998). 
LSND collaboration, Phys. Rev. Lett. 81, 1774 (1998). 



For a comprehensive review, see P. Bhattacharjee and G. Sigl, astro-ph/9811011, Physics Reports in press. 



FIG. 1. Feynman graph showing the computation of G v , where m„ = @u^= 

FIG. 2. The ratio m„/mjv (Eq. (23)) as a function of M2 (in units of Mf, and hence the notation Mf = 1), for Mp = 5 
and for various values of Ma- For visibility purpose, a few curves have been inflated by factors xlO 2 ' 3 ' 5,6 . 
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FIG. 3. The ratio m v /m N (Eq. (23)) as a function of M 2 (in units of M F , and hence the notation M F — 1), for M P — 50 
and for various values of Ma- For visibility purpose, a few curves have been inflated by factors xlO 2 ' 3 ' 5 ' 6 . 

FIG. 4. The ratio m v /m,N (Eq. (23)) as a function of M2 (in units of Mf, and hence the notation Mf = 1), for M p — 500 
and for various values of Ma- For visibility purpose, a few curves have been inflated by factors xlO 2 ' 3,5 ' 6 . 

FIG. 5. The ratio m^/mjv (Eq. (23)) as a function of M2 (in units of Mf, and hence the notation Mf = 1), for Mp — 5000 
and for various values of Mq. For visibility purpose, a few curves have been inflated by factors xlO 2 ' 3,5 ' 6 . 



FIG. 6. The ratio Ri = AI(G, P,Xi)/AI{G, P) for b = 0.035 



FIG. 7. The ratio i? 7 = AI(G, P, Xi)/AI(G, P) for 6 = 0.000095 

FIG. 8. The median mass m,2 as defined by Eq. (49b). Notice the correlation with m| — m 2 and m 2 — ml shown in the next 
two figures 



FIG. 9. m§ - mj as defined by Eq. (52a) for 6 = 0.000095 



FIG. 10. ml - ml as defined by Eq. (52b) for b = 0.000095 



FIG. 11. Diagram for m.34,43 



TABLE I 


. Particle content and quantum numbers of SU(3) C <g> SU(2) L i 


5C/(l)y ®SO(N f )®SU(2) lJR 


Standard Fermions 


9i = (3,2,1/6,^,1) 
II = (1,2,-1/2,^,1) 

uh = (3,1,2/3,JV>,1) 
d fl = (3,1,-1/3,^,1) 

c fl = (l, 1,-1, N f ,l) 


Right-handed z/'s 


Option 1: = (1,1,0, AT/, 2) 
Option 2: r? fl = (1,1,0,^,2); 

i/ fl = (l,l,0,l,2) 


Vector-like Fermions 


F£,h = (1,2,-1/2,1,1) 

MlL,fl = (1,1,-1,1,1) 
M2L,B = (1,1, 0,1,1) 


Scalars 


n a = (1,1,0,^,1) 

= (1,1,0,^,2) 
0= (1,2, 1/2, 1,1) 
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